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ABSTRACT 

We develop a framework based on energy kicks for the evolution of high- 
eccentricity long-period orbits with Jacobi constant close to 3 in the restricted 
circular planar three-body problem where the secondary and primary masses have 
mass ratio /i -C 1. We use this framework to explore mean- motion resonances 
between the test particle and the secondary mass. This approach leads to (i) 
a redefinition of resonance orders to reflect the importance of interactions at 
periapse; (ii) a pendulum-like equation describing the librations of resonance 
orbits; (iii) an analogy between these new fixed points and the Lagrangian points 
as well as between librations around the fixed points and the well known tadpole 
and horseshoe orbits; (iv) a condition a ~ yU^^/^ for the onset of chaos at large 
semimajor axis a; (v) the existence of a range /i < ~5 x 10~^ in secondary mass 
for which a test particle initially close to the secondary cannot escape from the 
system, at least in the planar problem; (vi) a simple explanation for the presence 
of asymmetric librations in exterior 1 : resonances. 



Introduction 



The three-body problem, or the dynamics of three masses due to their mutual gravita- 
tional influences, has a number of well-known special cases. One of these, the circular planar 
restricted case, requires that the primary and secondary bodies, mi and m2, follow circular 
orbits about their common center of mass and that the third body be a massless test particle 
moving in the massive bodies' orbit plane. These conditions simplify the three-body problem 
enough to produce an integral of the motion: the Jacobi constant Cj = —2{E — Qh) where E 
is the particle's energy^, h is the particle's angular momentum, and fl is the massive bodies' 
constant angular velocity. 



^We refer to the test particle as 'the particle' and to its energy per unit mass and angular momentum 
per unit mass as its 'energy' and 'angular momentum'. 
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Still, the circular planar restricted case has important applications to the dynamics of 
our solar system. Many of the orbits of major planets about the sun are nearly circular 
and are roughly confined to a plane; the same goes for many of the orbits of large moons 
about their planets. Common examples of applications for the circular planar restricted case 
include the effects of Jupiter on the asteroid belt; of Neptune on the Kuiper belt; of moons 
on planetary rings; and of giant planets on comets. 

This paper describes a study of this problem in the regime where m2 <^ mi, the test 
particle's eccentricity is large, and the Jacobi constant is greater than but close to 3 in 
the standard system of units where G = 1, the primary-secondary separation is 1, 1 = 
mi + 1712 — mi, and, therefore, Q = 1. Since values of Cj near 3 correspond to test particles 
on circular orbits close to the secondary, this special regime includes particles originally in 
circular orbits around a star close enough to a planet for the planet to perturb them into very 
eccentric orbits. Our interest in this regime arises from an intent to investigate the paths 
through which small particles are perturbed by a planet until they escape or are captured. 
This problem was studied by Ford et al. (2001) and Rasio & Ford (1996) via numerical 
simulations of three massive bodies in three dimensions. Due to this motivation we use 'star' 
and 'planet' to refer to the primary and secondary in the remainder of this paper. 

In §2 we derive to first order in /x = 7712/ {mi + m2) = m2 the energy kick received by 
a particle in a highly eccentric orbit with semimajor axis a ^ 1 at each periapse passage. 
We show that since the interaction is localized at periapse, this energy kick is essentially 
independent of a and depends only on the periapse distance and the azimuth difference 
between the planet and particle at periapse. In §3, §4, and §5 we use these energy kicks to 
find 'fixed' particle orbits and describe motion near them. These 'fixed' orbits are located at 
planet-particle mean-motion resonances. When observed stroboscopically at periapse only, 
they appear as fixed points just like the well-known Lagrangian points. We use both a 
continuous approximation and a discrete mapping in a derivation of the particle's motion 
around resonances, the resonance widths, and the libration periods. When these librations 
are observed stroboscopically, they likewise become analogies of the well-known tadpole and 
horseshoe orbits. In §6 we discuss types of chaos for large-eccentricity orbits, and in §7 we 
summarize and discuss our findings. 

2. Energy kick to first order in /j, 

Let AE the changes in the particle's energy between its consecutive apoapse passages. 
In our units, where the angular velocity of the star-planet is set to unity, the change in 
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angular momentum is also A£' ^. Therefore, it can be calculated by integrating the torque 
exerted on the particle: 

dV_ 

df 

where V is the gravitational potential produces by the planet and star and / is the particle's 
azimuth in inertial space. 



AE 



dt (1) 



To begin with, we estimate the energy kicks to first order in /x. Wc express AE as 
AE = fiAE^ + 0{i.P). To evaluate AE^ we calculate the torque assuming that the particle 
moves on a Kcplcrian trajectory around the star, with its focus at the center of mass. The 
effect of the deviation of the trajectory from that description on AE is of order ij? or higher. 

Since we are considering only the time elapsed between two consecutive apoapses, we 
choose coordinates such that the time t = when the particle is at periapse and the direction 
of periapse is along the positive x-axis. The planet and star are in uniform circular motion, 
so we can write V = V{9,r) where 9 is the angle between the planet and the particle and r 
is the particle's distance from the origin. This gives 



AE 

V is given explicitly by 



dV 



dt . (2) 



r 



y — 14>lanet + Ktar — , _ ^ T + , _ T \ (3) 

I '"planet | y ''star | 



to first order in /x, this gives 



1 f 1 cos^' 

^ ^ ^ ^ V + 1 - 2r cos ^) V2 ~ ^ J • ^> 



Let be the angle between the planet and the particle at periapse"^, so that 9 = (f)+t — f. 
Then the derivative with respect to 9 at fixed r which appears in Eq. (2) can be replaced by 
a derivative with respect to (j). To first order in /j, the particle trajectory r{t) can be assumed 



^Since the angular momentum is always perpendicular to the orbit plane, only one of its components 
is nonzero. We therefore treat the torque and the change in angular momentum as scalars equal to the 
components of the vector torque and the vector change in angular momentum which are perpendicular to 
the orbit plane. 

^Thus defined, cj) is the usual resonant argument measured at periapse only. 
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fixed and independent of 0, so we can move the derivative outside the integral of Eq. (2). 
Using the first order expression for V we get 



where the effective potential is given by 



(5) 



+ 1 - 2r cos 61 ^/r^ + 1 + 2r cos(i - /) 



+ 



cos 9 — cos{t — f) 



dt . (6) 



The integral is performed over one Keplerian orbit of the particle. 



In this expression for U^, the first term in the brackets is the 'direct' term; it represents 
the planet's contribution. The second term docs not contribute to AE'^; it keeps from 
diverging when a — > cxo and is obtained from the first term by substituting = vr. The third 
term is the 'indirect' term; it represents interactions with the star. AE'^ and its effective 
potential are functions of n, 0, and the particle trajectory shape, which determines r and 
/ as a function of t. Note that up to a constant, the effective potential is simply the 
time-integrated potential over the trajectory of the particle. 

When the apoapse distance a(l + e) is much larger than both 1 and the pcriapsc distance 
Tp = a(l — e), the perturbing effects of the star and planet on the particle near periapse 
dominate over perturbing effects on the particle elsewhere in its orbit. In this regime, the 
entire energy kick AE occurring between consecutive apoapse passages can be thought of 
as a discrete event associated with a particular periapse passage. In the limit as a diverges 
due to energy kicks but Cj remains constant, e — > 1, approaches a constant, and except 
near apoapse the entire trajectory approaches a parabola independent of a: a — > oo. If the 
particle is outside the planet's Hill sphere. 



lim Tr, 



lim 



2 \ a 



CI 



r = 2rp/ (1 + cos /) 

df _ (1 + cos/)^ 
dt 



(2r,)3/2 

t = (2rp)^/^ • ^ tan ^ + tan 



2/ 
2 



(7) 

(8) 
(9) 

(10) 



For particles that start close to the orbit of the planet, the periapse distance is therefore 
rp = 9/8. 
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Given this asymptotic form for the orbit, we can calculate the asymptotic forms of 
U^{a,e,(l)) — > U^{rp,(()) and A£'^(a,e, 0) — > A£'^(rp,0) in the large-a limit. For Cj = 3, 
the computed values of and its derivative AE^ as a function of are shown in Figures 1 
and 2. Near = 0, AE^ is dominated by the direct particle-planet interaction because the 
minimum planet-particle distance is much less than the star-particle distance. When = 0, 
AE^ = because of symmetry. When < but |0| -C 1, the planet lags the particle for 
most of the time the particle spends near per lapse, so AE^ < 0. Similarly, when > and 
101 < 1, AE' > 0. 

When 1101 — 7r| <^ 1, the indirect contribution AE^^^^ due to the star's reflex motion 
dominates because the star passes closer to the particle. From the 0-dependent part of star's 
contribution to the integral in Eq. (6), AE^^^ is a sinusoidal function of 0: 



where in the last step we use = 9/8 as an example in evaluating the coefficient. 

The integral in Eq. (11) seems to suggest that star-particle interactions over all intervals 
in / should contribute significantly to AE^^^. However, as Eq. (10) shows, \t\ increases much 
faster than |/| as |/| approaches vr. As a result, oscillations in cos(t — /) kill contributions to 
the integral at |/| near tt, and the star-particle interaction is important only near periapse. 

We can also get the total contribution of terms that are second order or above in this 
is just the difference between values of AE found by numerical integration of the equations 
of motion and values of /lAE^ given by Eq. (6) (see Figure 1). 



Resonances occur when the particle completes p orbits in exactly the time needed for 
the planet to complete p + q orbits for some integers p, q. This situation is known as a 
p : p + q resonance. In the standard treatment of these resonances, both orbits in question 
are usually nearly circular and a significant interaction occurs every time the bodies are 
at conjunction, i.e. whenever their azimuths coincide. This happens once every resonant 
cycle if g = 1, so g = 1 resonances arc usually termed 'first order' resonances. During a 
conjunction between a test particle and planet in orbits with low eccentricity e -C 1, the 
torque exerted on the particle while the particle precedes the planet almost cancels the 
torque exerted while the particle lags the planet; the residual is of order e. When q > 1, q 
conjunctions occur during each resonant cycle. Because they occur in different positions in 
inertial space, their effects tend to cancel each other, leaving a residual torque of order e*. 




(11) 



3. 



First-order resonances 
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Fig. 1. — First-order energy kick AE^ in the large-a regime (solid line) with Cj — 3. The 
dotted line is the energy kick AE/pi calculated for a = 10~^, Cj — 3 parabohc orbit with 
all higher order terms included. For this /i, the first-order term clearly dominates for all 
values of 0; higher-order effects in are visible only near = —0.12. The dashed line is the 
planet's direct contribution to AE^; the dash-dotted line is the indirect contribution to AE^ 
from the star's reflex motion. 

Since the interaction strength decreases exponentially with increasing g as e*, resonances in 
the standard treatment are usually classified by q value. Accordingly, a p : p + q resonance 
is called a 'gth-order' resonance regardless of the value of p. 

However, the high eccentricities of orbits in the large-a regime discussed here make the 
standard definition of resonance order meaningless. Since e — > 1, resonances of different 
order under the standard definition have comparable significance because ~ 1. Also, 
encounters at periapse are physically far more important than conjunctions at other points in 
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Fig. 2. — Potential in the large-a regime with Cj — 3. 

the particle's orbit. We therefore redefine the 'orders of resonance' to focus on interactions at 
periapse. If the planet completes an integer number of orbits in the time it takes the particle 
to orbit exactly once, then we say the particle is in a 'first-order' resonance. In general, 
we say the particle is in an 'pth-order' resonance if the planet completes an integer number 
of orbits in the time it takes the particle to orbit p times: then there are p interactions 
within one resonant cycle. In terms of the standard resonance treatment, we say Sk p : p + q 
resonance in the large-a regime is 'pth-order' regardless of the value of q. In both the large- 
and small-eccentricity cases, the order of the resonance is given by the number of significant 
interactions within a single resonant cycle. 

In the following we show that this revised definition does indeed make sense. We calcu- 
late the widths of resonances of various orders in the large-a limit and show that with this 



-8- 



new definition, their widths decay exponentially with the order of the resonance. We discuss 
in detail the first-order or 1 : A^" resonances. 

According to this definition of resonance orders, should be constant in time if we 
consider a particle exactly at a first-order resonance of semimajor axis Ores and if we ignore the 
effects of energy kicks. A particle close to resonance with, say, semimajor axis a — Ores + Aa 
should drift in over time at a constant rate, again ignoring energy kicks. The amount of 
drift per orbit of the test particle is just the difference between its orbital period 2na^^'^ and 
the resonant one 27raUs- We can express this drift as 

^ = 3W/2Aa . (12) 

an 

Here the dn acts as a place-holder indicating that the d(f) is associated with a single particle 
orbit. The differential is justified assuming Aa <^ a~^/^ so that the change in per particle 
orbit is much less than tt. We refer to this differential form as the continuous approximation. 

Energy kicks cause the semimajor axis to evolve in time. To first order in /i we have 

d(Aa) ^dE ^ ^ dU^ 
dn dn d<p 

To justify the differentials here we require that /x be small enough for the change in Aa due 
to a single kick to be much less than the typical Aa. We differentiate Eq. (12) and substitute 
Eq. (13) to get 

« 5/2 dU^ nA^ 
This shows that (f) simply evolves as a particle moving in the potential U^{4>)- 



3.1. Generalized Lagrangian points 

Since has four extrema at the four zeroes of AE^, there are four fixed points. These 
fixed points are simply periodic orbits of the particle. Two of the fixed points are unstable 
since they correspond to maxima of the potential 9 — 0,7r. The other two are stable since 
they correspond to potential minima at ^ = ±1.21. The existence of two extrema at ^ = 0, tt 
is guaranteed by symmetry arguments. The two additional extrema at (f) — ±1.21 occur 
where the energy kicks from the planet and star cancel each other exactly. These extrema 
therefore appear only when the indirect term — or, equivalently, the star's motion — is taken 
into account. 

This discussion suggests an analogy between the five well known Lagrangian points and 
the new fixed points. The two stable points correspond to the stable Lagrangian points 
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L4 and L5, which also appear only when the motion of the star, i.e. the indirect term, is 
taken into account. The unstable fixed point at = tt is the analogue of L3; the one at 

= corresponds to Li and L2, which merge in this generalization. For a given resonance 

1 : N, N = (T'l'^, wc therefore denote the fixed points by L^, L^, ^ and . The 
positions of these new fixed points in comparison to their standard Lagrangian counterparts 
is summarized in Table 1. 





Lagrangian points 


generalized Lagrangian points 


resonant index 


(iV = l) 


iV = 2 


iV = 3 


iV = 4 


large a 


semimajor axis 


1 


22/3 


32/3 


42/3 


a = iV2/3 


physical meaning 


particle is stationary 


particle moves 


on periodic orbit 


of fixed points 


in rotating frame 




in rotating frame 




definition of 


azimuth of particle 


azimuth of particle in rotating frame 


angular variable 


in the rotating frame 




when it is at periapse 






01 = 0, 02 = 


single point with 0i2 = 


= 


is 


03 = TT 




03 


= TT 




-L4 & -^^5 (05 = -04) 


</,4 = 7r/3 ~ 1.04 


0| = 1.196 


0| = 1.196 


(pi = 1.198 


04 = 1.21 


min. tadpole period 


^M-'/' ^ 2.42/.-V2 


4.4//-V2 


5.1//-1/2 


5.5^-1/2 




Aomax tadpole 


^/8/3n^/'^ ~ 1.63//V2 


1.4//V2 


I.6//I/2 


I.8//I/2 


0.78a3/V'/' 


Aamax horseshoe 


2(3)1/6^1/3 ^ 2.40/xV3 


4.6AiV2 


4.7/xV2 


5.0/xV2 


l.Sa^/V^/' 



Table 1: Comparison of properties of Lagrangian points Li_5 and orbits around them with 
those of their generalized versions ^^^^5. All quantities are given to lowest order in /x. In 
particular, expressions for the iV = 2, 3, 4 resonances were calculated using a potential com- 
puted to first order in at a = N'^^^ rather than in the large-a limit. The numerical values 
for the generalized Lagrangian points and orbits are given for Cj = 3. 



3.2. Generalized tadpoles 



The analogy is more obvious when motion around the fixed points is investigated. Small 
amplitude motion around the stable fixed point and can be approximated by expand- 
ing around its minimum. This results in a harmonic oscillator equation: 



d-n? 



(0 - 0r 



(15) 
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The small amplitude libration period around either or is therefore given by 

K = ^ii^ =(^^ V «-'/V-^/' = 0.79a-^/V-^/^ (16) 

where in the last step we use Vp — 9/8 in the large-a limit to get (PU^/dcf)'^ ~ 3.3 at = 04,5. 
Note that K gives the number of periapse crossing per libration period. In our units, where 
27r is the period of the massive bodies, the libration period is then 27ra^/'^K. 

Since Eq. (14) describes motion under the influence of a fixed potential, we can write 
down the conservation of energy equation by multiplying Eq. (14) by ^ and integrating with 
respect to n: 

-(^) + 67ra5/Vt^^ = constant . (17) 

The constant of integration is the 'energy' associated with the movement of the orbit in 
and a. Since the potential is finite, it can only support a finite particle 'speed' in libration 
around or . The 'speed' is directly related to the deviation of the semimajor axis from 
the resonance via Eq. (12), so the maximal width of these librations in a is given by 

/ 4 \ V2 

Aa_ = (^i- j /.V2^3/4[^i(^) _ c/i(0^)]i/2 ^ o.78a^/V^/^ ■ (18) 

These librations around the fixed points or are analogues of the well-known 
tadpole orbits. Note that the maximal widths of both the standard and generalized tadpole 
orbits scale as //^/^ (see Table 1). The similarity is more apparent if we treat the (a, 0) pa- 
rameters, which describe the orbit of the particle, as polar coordinates as shown in Figure 3. 
Seen in this way, (a, 0) are analogous but not identical the polar coordinates of the particle 
in the rotating frame: a is the semimajor axis, not the radius, and is the azimuth of the 
test particle in the rotating frame only at periapse passage. Then the fundamental difference 
between the (a, 0) plane and the rotating frame is that while generalized Lagrangian points 
and the motion around them exist in a surface of section made up of discrete points repre- 
senting periapse passages, the standard rotating frame with the standard Lagrangian points 
is made up of continuous trajectories. Therefore, while the standard Lagrangian points are 
fixed points in the rotating frame, the generalized points represent periodic orbits in that 
frame. 

The generahzed tadpoles are equivalent to 'asymmetric librations' — trajectories whose 
resonant argument librates about a value other than or tt. In this context, and 
correspond to 'asymmetric periodic orbits' whose resonant argument is constant but not 
equal to or tt. Our discussion above gives a simple physical argument for the existence of 



-2 





a cos (f) 



2 



Fig. 3. — The Lagrangian point analogues Lf for N — 1,2,3,4. The diagonal lines trace 
the azimuths of Lf^^ — that is, the (f) values of the minima in U^. Sample 'horseshoes' and 
'tadpoles' calculated via numerical integration with = 2.5 x 10~^ and Cj — ?> are also 
shown. 

asymmetric librations in all stable 1 : N exterior resonances. Again, note that the existence 
of these asymmetric librations and asymmetric periodic orbits follows from analysis of 
only when both the direct and indirect terms are accounted for. 
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3.3. Generalized horseshoes 

As the energy of the particle moving under the potential increases beyond that of 
the maximal tadpole orbit, it overcomes the lower potential barrier at = tt. As long as its 
energy is still below the higher barrier aX <p = 0, the particle will librate around both the 

and points, avoiding only a narrow range in around = 0. These trajectories are 
the generahzed horseshoe orbits. Using the same method as we used for the tadpoles, we 
calculate their widths in the continuous approximation to be 

Aa^ax = ^i'/'a'/'[U\0) - U\cl>,)]'/' = l.Sa'/Y^' . (19) 

The width of the maximal standard horseshoe does not follow this ^^/^ pattern, since the 
standard horseshoe case differs qualitatively from its generalized version. For the standard 
horseshoe, the angular momentum change is concentrated near the horseshoe's two ends. The 
close approach of the particle to the planet there increases the strength of the interaction 
beyond As a result, the width of the horseshoe scales as //^/^ rather than //^/^. For a 
generahzed horseshoe, the hbrating particle never gets closer to the planet than — 1. 

In Figures 4 and 5 we show libration around L45 and L\^^. In order to focus on the 
motion close to these points we plot a and as Cartesian rather than polar coordinates. In 
these plots, the librations in the surfaces of section appear to be 'warped' when compared to 
the continuous approximations calculated using the pendulum-like Eq. (14). This 'warping' 
is due to the discrete nature of the motion in the surfaces of section. As a trajectory moves 
from Aa = toward larger positive Aa values, for example, the energy kicks stay positive 
and Aa should keep increasing until the trajectory reaches a </> value corresponding to a zero 
in the AE'^ vs. curve. Within the continuous approximation, we expect the trajectory to 
begin moving back toward Aa = at exactly this (p because AE'^ changes sign. A discrete 
trajectory will 'overshoot' the nominal (j) where AE^ — since a positive energy kick will 
carry the trajectory past this before the first negative kick is apphed. As a result, the 
libration trajectories in the surfaces of section tend to become warped in the direction in 
which orbits move when librating. A quantitative discussion of this feature is given in the 
next section. 



4. The eccentric mapping 

The 'warping' noted above suggests that the discrete nature of the surface of section is 
essential to understanding some feature of the motion in the (a, 0) plane. To study this, we 
define a mapping from the (a, 0) plane to itself. Beginning at an arbitrary point, this mapping 
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2.48 2.52 2.56 
a 




Fig. 4. — a vs. (f) plot for N = 4, fi = 10~^, Cj = 3. We use [0, 27r] as the range in to 
show tiie trajectories more clearly. The left-hand plot contains trajectories computed under 
the continuous approximation. The middle plot contains a surface of section computed via 
numerical integration of the orbits. The right-hand plot contains trajectories computed via 
the eccentric mapping discussed in §4. The same initial conditions were used for the tra- 
jectories in all three plots. The continuous approximation plot lacks the chaotic behavior 
evident in the numerical integration and eccentric mapping plots. Trajectories in the map- 
ping plot differ from the numerical integration plot mostly because they were calculated with 
U^, the potential in the large-a limit. Note that the separatrix trajectory in the middle plot 
is chaotic but on a scale too small to see in this figure (see Figure 9). 



produces an infinite series {a\ 0*) of points visited by the test particle in the (a, (p) plane. 
Except perhaps in the few lowest- resonances, we can build an excellent approximation to 
the correct mapping by applying the first-order kicks in the large-a limit: 

--^ ^ --^ + liAE\(P') (20) 



0*+i=0^ + 27r(a*+i)3/2 



(21) 
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Fig. 5. — Same as Figure 4 but for N = 10. For this larger A^, the resonances are wider. So 
the resonance overlap is more severe for the outer edges of the = 10 resonance than for 
the N — A one. This causes the destruction of all horseshoe orbits and the distortion of the 
tadpoles relative to those computed in the continuous approximation. 



where the new value of is calculated modulo 27r, i.e. brought back into the interval (— tt, tt) 
by adding an integer multiple of 2tt. Note that the a-value used to find (ff^^ itself has 
index i + 1; physically, this corresponds to the large-a limit assumption that each energy 
kick is a discrete event associated with a given periapse passage. Applying this mapping 
for several initial values in the (a, (p) plane results in the right panels of Figures 4 and 5. 
The close resemblance between trajectories generated with the mapping and with numerical 
orbit integrations demonstrate this mapping's accuracy. 

It turns out that the warping of the small amplitude tadpoles can be understood 
completely in terms of the mapping. Close to the fixed points and L^, we define 
Aa' — — ares and A0* = 0' — Ores so that the mapping becomes 

(22) 



Aa^+^ = Aa* - 2a',^^nA(f)' 
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(23) 



Since these are linear recursive equations, they can be solved analytically by standard tech- 
niques. We seek a solution of the form (Aa*, A0') = (^4, $)q;\ (The reader should not 
confuse indices with exponents — superscripts above Greek letters other than are expo- 
nents.) Substituting in the recursive equations, and seeking a non-trivial solution we obtain 



a 







(24) 



P4,5 



Note that the dimensionless parameter in this equation is simply {27i/K)^ where K is the 
number of periapse passages per libration in the continuous approximation as given by 
Eq. (16). If we denote the solutions as ai and a2, it is clear from the above equation 
that their product aia2 equals 1. Since we are interested in potential minima, > 0. 

For K > 71, the two roots are complex conjugates and each has unity norm. The fixed 
point is therefore an elliptical point in the discrete mapping as well as in the continuous 
approximation. The two values of a are given by 



ai^2 = 1 - '^{T^/Kf ± 2^1 - (7r/X)2(7r/ir)i 
The number of periapse passages per libration is given by 



^map 



27r 



arg(Q;) 



27r 



arctan 



'2 ^l-{n/Kf{n/K ) 
1 - 2(7r/X)2 



(25) 



(26) 



Clearly, as X — > oo, i^T^^^P/X — > 1. This is expected since the continuous approximation is 
justified in this limit. Using the two values of ct, we can find the eigenvectors: 



(A$) 



, {T:/Kf ± Vl - {n/Ky{n/Ky, 



(27) 



Since the eigenvectors determine the axes of the ellipses representing small librations about 
the fixed points, the similar shapes and orientations of the smallest librations in the middle 
and right-hand panels in each of Figures 4 and 5 confirm the eccentric mapping's accuracy. 

For a = lO^/^, fi = 10~^, the continuous approximation gives 11.6 orbits per tadpole 
libration (Eq. 16). The eccentric mapping gives 11.4 orbits (Eq. (26)). This is close to the 
the 10.7 orbits per libration observed for very small librations about the fixed points'^. The 



"'The largest tadpole libration shown in Figure 5 breaks into 14 islands. This is an example of the 
Poincare-Birkhoff fixed-point theorem. It indicates that this tadpole's libration period is 14 orbits. This 
lengthening of the period is expected as the trajectory grows toward the separatrix passing through (p = it. 
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negative power of a in Eq. 16 implies that as a increases, the number of periapse passages 
per tadpole libration period will decrease and the trajectory shapes will become increasingly 
warped. 

In fact, when a grows so large that K falls below tt, the tadpoles are destroyed. For 
X < TT, the roots of Eq. (24) are real and distinct; therefore one of them is larger than unity. 
Then the fixed point is not stable despite being at a potential minimum. Our quantity K is 
closely related to the residue R discussed by Greene (1979): R—1 — {tx/KY. 

The warping of the tadpoles which, at its extreme, leads to destruction of the resonances 
is absent in the continuous approximation. However, it can be understood as perturbations 
from nearby resonances. Interactions between neighboring first-order resonances should be- 
come large enough to destroy these resonances when the resonances begin to overlap. Eq. (19) 
implies that as a increases, the resonances widen in a while the distance between them de- 
creases. Then we can find a condition on and a for resonance overlap by dividing half the 
distance between consecutive first order resonances by the width Aamax of each resonance 
as given by equation (19). In the large-a limit, the distance between resonances is given by 
ia"-*^/^ so we obtain 



3 



resonance separation / tt ^^ V2 



2Aa_ -^^2) a-^V^/^r(O)-C/^(04)]-^/^ = O.18a-VV-V^ . (28) 

This is proportional to K^: in the large-a limit with Vp = 9/8, the right-hand side is 0.28K'^ 
and first-order resonances overlap when < 1.9. In this case, therefore, first-order reso- 
nances are destroyed before they formally overlap. With /j, = 10~^, first-order resonances 
are destroyed above a ~ 4.02. Then if we neglect higher-order effects in fj, and differences 
between the potentials for a ~ 4 and the large-a limit, we expect just eight stable first-order 
resonances. 



5. Higher-order resonances 

As defined in §2, higher-order resonances are the p : p + q resonances with p > 1. These 
resonances are located at ares = (iV/p)^/^ where N — p + q is a,n integer relatively prime to 
p. In analogy to our treatment of first order resonances, we note that if we neglect energy 
kicks, a particle exactly at resonance should move in (p by 27r during each resonant cycle and 
by 2'nq/p between consecutive periapse passages. The stationary points of this resonance 
should therefore occur at regular intervals of 27: /p in (j). 

To study motion near but not at resonance, we include energy kicks. For a particle close 
to resonance, we can follow its trajectory by treating each resonant cycle as p applications 
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of the eccentric mapping, one for each periapse passage in the cycle 



i=0 



0=9i^+27rj/p 



37raJe4^ 1 - 67rar/s^/i 



(29) 



As before, Aa^ — — a^es and A(j>' = (fP — (p^ where (p^ corresponds to the nearest fixed 
point in the resonance. The condition under which the linearization in is valid is now 
Aa^ <ti a^/^/p^ instead of Aa^ <^ a^l'^ because the number of energy kicks per resonant cycle 
is "p instead of 1 and because the scale in over which the potential changes is now vr/p 
instead of tt. The condition under which linearization in [i is valid also changes because the 
largest term linear in \i that appears in the mapping matrix is aUsl^^^r- Though /i itself 

is small, cross-terms of order /x^ and higher are now important unless (his 1^^^^ <^ 1- This 
stronger condition is equivalent to X ^ 1, so the higher-order resonance treatment does 
not offer any simplifying advantages over the eccentric mapping discussed in §4 unless K is 
large. 

Since changes very little between consecutive periapse passages in this K 1 
regime, we can use a variant of the continuous approximation where we neglect the effects of 
drift in within a single resonant cycle. Then we can treat the particle's motion in terms of 
the net energy kick over an entire resonant cycle rather than a single particle orbit. The net 
energy kick is just the sum of p energy kicks spaced 27r/p apart in 0, so the particle appears 
to move in the potential 

p-i 

Ul = Y,U\<t>-27:k/p) . (30) 

fe=0 

Note that effects of the star's reflex motion do not contribute to C/^ if p > 1: the indirect term 
in is exactly sinusoidal and the sum of p identical sine curves spaced 271 /p apart in phase 
is 0, so f/p iiid = 0. Since the part of due to the planet's direct contribution has just one 
maximum and one minimum at = 0, tt respectively. Up has p identical maxima and minima 
(see Figure 6. Then a trajectory librating in one of the minima of should appear as a series 
of 'islands' spaced evenly in in the (a, 0) plane. As a result, no asymmetric librations are 
possible in higher-order resonances. Our result that, among exterior resonances, only 1 : N 
resonances show asymmetric librations is consistent with work done by Prangakis (1973). 
He analyzed expressions for the time-averaged direct and indirect terms of the disturbing 
function to find that asymmetric librations can exist only in p '■ P + q resonances where 
p = ±1. 

Because the p energy kicks received by the particle during each resonant cycle are spaced 
evenly by 27r/p in 0, we expect that the kicks will partially cancel over each resonant cycle 
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Fig. 6. — U^Q, or Up for a lOth-order resonance: summed over 10 consecutive periapse 
passages spaced evenly in 4>. As before, Cj — 3. 

and that this cancellation will improve exponentially as p increases. We therefore expect the 
amplitude of Up to decrease exponentially with p. As Figure 7 shows, this exponential decay 
is observed numerically: a best-fit line in log-log space gives amplitude oc 1.20~^. 

6. Chaos in the Icirge-a Hmit 

We discuss just a few of the types and regions of chaos which arise when a is large. 
We first discuss 'global' chaos, which consists of chaotic regions that span a few resonance 
widths or more. We then give a few examples of 'local' chaos — chaos confined to regions 
within a single resonance — and compare the structure seen in trajectories in the (a, 0) plane 
on local and global scales. 

On large scales in a, chaotic regions arise where there is overlap between neighboring 
resonances or instability due to a small winding number as discussed in §4. In regions of the 
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Fig. 7. — Amplitudes of Up plotted on a log scale as a function of p. Used Cj — 3. Best-fit 
line is log^g [amplitude] = — 0.078p + 0.21, or amplitude oc 1.2~^. 

(a, (f)) plane where first-order and/or higher-order resonances overlap even partially, we expect 
to see contiguous 'globally' chaotic regions which span large ranges in a. Any remaining 
stable regions within resonances will appear as 'islands' of stable librations. Particles can 
undergo large changes in a only if they move in these chaotic regions, so such regions provide 
the only channels through which initially bound particles can escape from the star-planet 
system. 

If a is large enough, K falls below tt and these 'islands' disappear as discussed in §4. 
For a given value of /i, we see from Eq. 16 that this occurs when 



where, again, the numerical example corresponds to = 9/8. The condition K < tt for 
resonance destruction in higher-order resonances does not follow simply from Eq. 28: effects 
of order //^ or higher may make the winding number expressions for higher-order resonances 
differ from the first-order resonance case in Eq. 16. However, numerical experiments suggest 
that the a at which higher-order resonances become unstable is comparable to but less than 




(31) 
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that given by 31. 

If fjL is small enough, there should be regions in a where the resonances do not overlap. 
In these regions we expect to see stable trajectories which circulate around the resonances 
instead of librating in them. We find numerically that stable circulating trajectories exist 
for fjL values up to at least ji — b x 10~^; an example is shown in Figure 8. Greene (1979) 
suggests that as increases, the last stable circulating trajectory should have semimajor 
axis a such that a^/^ is the golden ratio (1 + v^)/2. Our situation differs qualitatively 
from Greene's in that our potential depends on its linear coordinate, the semimajor axis, 
while Greene's potential, which is given by the standard map, is independent of its linear 
coordinate r. Specifically, when a is not much larger than 1, — 1 ^ 1; this leads to a 
larger maximum energy kick and potential well depth than is expected for = 9/8, so the 
resonances are wider and more prone to overlap for a given a close to 1 than we would expect 
in the large- a limit. However, as a increases the resonance spacing decreases as discussed 
in §4. These competing effects suggest that the last stable circulating trajectories — those 
which, in a sense, are 'farthest' from any resonances — should lie neither near a = 1 nor 
at a 3> 1. Also, effects of order /i^ and higher which are present in our situation have no 
analogue in Greene's analysis of the standard map. So it is unsurprising that the last stable 
circulating trajectories which we found numerically have a^/"^ unrelated to (1 + \/h)/2. 

Continuity and uniqueness imply that in a system with two degrees of freedom, stable 
trajectories in a two-dimensional surface of section cannot be crossed. In the planar restricted 
three-body problem, therefore, stable circulating trajectories divide the (a, 0) plane into 
separated regions in a. This implies that for any //, < 5 x 10^^, chaotic and regular trajectories 
which start close enough to the planet are confined to a set range in a. Then the particles 
associated with these trajectories can never escape from the star-planet system. 

This bounding of chaotic regions by stable trajectories also leads to confinement of chaos 
on very small scales in a. Regions of small-scale 'local" chaos arise from unstable fixed points 
which must be saddle points due to the area-preserving nature of the eccentric mapping; the 
separatrices associated with the saddle points are chaotic. If stable continuous trajectories 
exist near a saddle point, they act as boundaries to the chaotic separatrix. Prominent exam- 
ples of these separatrices include those dividing tadpole and horseshoe analogue trajectories 
within individual first-order resonances. These regions are bounded by the largest stable 
tadpole and smallest stable horseshoe, so their maximum range in a is at most the resonance 
width. One of these is shown in Figure 9. 

The existence of similar separatrices on all scales smaller than a single resonance width 
follows from the Poincarc-Birkhoff theorem, which states that for small enough yU, a trajectory 
with rational winding number K is associated with equal numbers of alternating stable and 
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Fig. 8. — A stable circulating trajectory in (a, 0) space calculated via numerical integration 
with = 5 X 10"^ and Cj = 3. 



unstable fixed points. According to the KAM theorem, some continuous trajectories — that 
is, trajectories with irrational K — should also be stable as long as is small enough. If 
a trajectory with rational K is bounded on either side by stable continuous trajectories 
with irrational then the chaos associated with the unstable fixed points is confined to 
the region bounded by the continuous trajectories. As for the stable fixed points, they are 
associated with their own librating trajectories; the tadpole analogue made up of islands 
shown in Figure 5 gives an example of such librations. We expect some librations like these 
to have rational winding numbers and, therefore, their own sets of unstable fixed points 
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and confined chaos on an even smaller scale. In principle, this argument can be applied 
repeatedly within a single resonance to unearth similar chaotic regions on scales as small as 
desired. 

We can treat the entire (a, 0) plane as an extension of this self-similarity to the largest 
possible scales. If we plot (a, 0) as polar coordinates, a p '■ P + Q resonance trajectory 
appears to 'wind' around the point a = with rational winding number p/ (p + q). Also, the 
corresponding resonance is associated with p stable and p unstable fixed points when p > 1 
and 2 stable and 2 unstable fixed points when p = 1. This provides a striking visual analogy 
to the librations seen within a single resonance. 

7. Discussion and conclusions 

Using simple physical reasoning instead of explicit analysis of terms in the disturbing 
function, we have developed a framework for studying particles perturbed into exterior high- 
ccccntricity orbits in the circular planar restricted three-body problem. We have found that, 
to first order in /i, these orbits move in (a, (j)) phase space according to a potential with max- 
ima at = 0, TT separated by symmetrical minima. In the special case of resonance orbits, 
movement in this potential translates into behavior governed by a modified pendulum equa- 
tion. Previous pendulum-analogue analyses of this problem have usually been formulated 
via the disturbing function and the continuous resonant argument (Winter & Murray 1997a; 
Dermott & Murray 1983). 

Our analysis, specifically that of mapping, is most similar to that of Malyshkin & 
Tremaine (1999). They consider the evolution of high-eccentricity comet-like orbits in the 
low-inclination circular restricted three-body problem by integrating numerically to find the 
energy kick as a function of the resonant angle at periapse and then using this energy kick 
to create a mapping which takes one periapse passage to the next. However, Malyshkin & 
Tremaine (1999) are interested in particle orbits that cross the orbit of the secondary, so the 
form of their energy kick is qualitatively different from ours. In particular, while a small 
nonzero orbital inclination would barely affect our energy kick function, it could drastically 
change the shape of the overall energy kick function in the case of planet-crossing orbits. 
Partly because of this, they do not discuss their energy kick it in terms of a potential. Also, 
they focus on the chaotic diffusion of the particle toward escape or capture rather than on 
motion in resonances. 

For 1 : N resonance orbits — that is, those we call first-order resonance orbits — the shape 
of the potential generates analogues of the Lagrangian points for > 1 resonances. The 



-23- 



potential similarly leads to two kinds of libration analogous to the horseshoe and tadpole 
orbits seen in a 1 : 1 resonance, p : resonances — that is, those we call higher-order 
resonances — show only one kind of libration: when the winding number is large, the sum 
relating the higher-order resonance potentials to the first-order resonance potential eliminates 
the indirect term responsible for the tadpole analogues. 

Several authors discuss the tadpole analogues' presence or absence in mean-motion res- 
onances in general under the name 'asymmetric librations'; Nesvorny & Roig (2001) are the 
only others we know of to refer to the 1 : resonance librations as 'tadpoles' and 'horse- 
shoes', though they do not elaborate on this analogy. Some authors have used analytical 
studies of the Hamiltonian and the disturbing function to set conditions for the existence 
of asymmetric resonances (Bruno 1994; Frangakis 1973; Message 1970). In particular, Pran- 
gakis (1973) analyzed the time-averaged direct and indirect parts of the disturbing function 
to deduce that only what we call first-order resonances should show asymmetric librations. 
Bruno (1994) also found analytically that asymmetric librations only exist in what we call 
first-order exterior resonances. We confirm this and provide a simple physical explanation. 

Others have used numerical methods to confirm the existence of asymmetric librations 
for particular 1 : N resonances and ranges in eccentricity (see, for example. Winter & Murray 
1997b; Beauge 1994; Message & Taylor 1978; Frangakis 1973; Message 1958). Some of these 
also compare their numerical results to expressions for the Hamiltonian correct to first or 
second order in eccentricity. Although the agreement is generally good for what we call 
first-order resonances, the Hamiltonian expressions for what we call higher-order resonances 
tend to predict spurious asymmetric librations. We believe these are due to extra extrema 
introduced into the potential by the cos (2 - resonant angle) term in the Hamiltonian. In some 
of the more recent studies involving asymmetric hbrations (Nesvorny & Roig 2001; Malhotra 
1996) the discussion is framed in terms of the dynamics of the classical Kuiper Belt and so is 
confined mostly to what we call low-A^ first-order and low-p higher-order resonances in the 
low- to moderate-eccentricity regime. 

We find a limit on a for stable first-order resonances. Overlap between the resonances 
creates chaotic regions of (a, 0) phase space; for semimajor axes larger than some a oc IJ>~^^^, 
the resonance centers are overlapped and no stable librations are possible. This is the high- 
eccentricity analogue of the well-known chaotic criterion |a — 1| ~ /x^/^ found by Wisdom 
(1980) for the circular planar restricted three-body problem in the low eccentricity case. 
We use the Chirikoff criterion for resonance overlap to estimate the location of the onset of 
chaos. For sufficiently narrow resonances, or small enough /i, there exist regions in (a, 4>) 
space which lie outside all of the resonances but which are not chaotic. In the planar problem 
we consider, particles interior to the circulating trajectories in these regions are never able 
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to escape from the star-planet system. 

The basic framework for the behavior of high-eccentricity orbits and the properties of 
chaotic regions in (a, 0) space can be apphed to the orbital evolution of small bodies in the 
solar system. Objects in the Kuiper Belt, for example, are believed to have arrived there 
via interactions with Neptune (see, for example, Malhotra et al. 2000); we can apply this 
framework to study their trajectories. Many of these objects are known to be in resonances 
(see Chiang et al. 2003, for a recent compilation). The mass ratio between Neptune and the 
Sun is = 4.4 x 10^^. Since this is above the critical fi ^ 5 x 10^^, Kuiper Belt objects 
with Cj = 3 are either librating around a resonance or moving chaotically. The latter could, 
in principle, be ejected as there is no stable circulation for that value of n = /j.^ and Cj = 3. 
However, the known Kuiper Belt objects span a range in Cj roughly 2.6 < Cj < 3.2. In the 
planar problem, = /xjv and for example Cj — 3.1, stable circulations exist, and protect 
some of these objects from escape. To study the ultimate fate of such Kuiper Belt objects, 
the effect of inclination must be understood. 

Similarly, we might apply this framework to the scattering of small planetesimals by 
giant protoplanets and could provide insight to numerical integrations such as those of Rasio 
& Ford (1996); Ford et al. (2001). Studies hke this require an investigation of the way 
in which the energy kicks move orbits through the 'global chaos' region surrounding the 
resonances in the (a, 0) plane. Although the antisymmetry of AE^{(f)) about = suggests 
that the orbits should random walk through phase space, effects of nearby resonances (e.g. 
Malyshkin & Tremaine 1999) and terms of higher order in /x become important on timescales 
long enough for escape to become possible. The importance of second order effects may be 
understood as follows. Since the amount of extra energy needed to escape is 1 / (Xjnit ^ 1 
and the energy kick per orbit is ~ //, we expect that the average number of kicks needed to 
escape is ~ Note that unlike the first order kicks, the 0{iJ?) energy kicks do not average 
to over the interval (— tt, tt] in 0. This is also apparent from Figure 1. Therefore, with //~^ 
kicks, the sum of 0(/i^) effects produced by the energy kicks will be of order unity — that is, 
of size comparable to the total first order effect. 

MP is supported by an NSF Graduate Research Fellowship and RS holds a Sherma 
Fairchild senior research fellowship. We thank Peter Goldreich for useful discussions. 
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Fig. 9. — A single chaotic trajectory corresponding to the separatrix dividing 'tadpole' and 
'horseshoe' hbrations in the N — A resonance when fj, — 10~^ and Cj — 3. This trajectory 
was computed via numerical integration with the same initial conditions as were used to 
produce the separatrix trajectory in the middle panel of Figure 4. It is confined in the (a, (f)) 
plane by stable librations similar to the smallest horseshoe and largest tadpoles shown in 
Figure 4. Note the empty spots in the outer reaches of the chaotic trajectory; these 'avoided' 
areas correspond to islands of stable librations around stable fixed points in trajectories with 
rational winding number. 



